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I. INTRODUCTION 



Quasi-exactly solvable (Q.E.S.) Schrodinger equations are of great current interest. Indeed they have the charac- 
teristic property that analytic expressions for a finite number of energy levels and their corresponding eigenfunctions 
can be obtained. These Q.E.S. equations have been first investigated in and since then they have been studied 
extensively (for a detailed review see f.i. |^). Actually one can find two different points of view for analyzing these 
equations. 

First, the number of analytic solutions is related jsj to the dimension of the irreducible representations of the 
Lie algebra sl{2,R). More specifically the corresponding Schrodinger Hamiltonians can be written, up to a change 
of variables and a change of eigenfunctions, as linear and quadratic combinations of the sl(2, R) generators, these 
preserving the finite-dimensional space of the solutions. 

Second, the number of analytic solutions is fixed at the start (essentially two Q or three ||]). Except in the case 
of two solutions the corresponding Schrodinger Hamiltonians cannot, in general, be related to any Lie algebra. 

Each of these two points of view has its own advantages and disadvantages. The Lie algebraic approach is elegant, 
straightforward and the associated potentials can be treated in a global way. However it concerns only a few interac- 
tions which have been listed in ^ . The second approach deals with a bigger and not yet known number of potentials 
but no general method to study them is available. Until now, the only method found in the literature to construct 
these Q.E.S. potentials is the one developed in Q| and It is based on a recursive construction of a finite number 
of supersymmetric partners . It will be recalled in Section [V with more details but let us already mention that the 



Kuliy-Tkachuk method is too constraining and thus misses some Q.E.S. interactions. Some of them are given in the 
following. 

The purpose of this paper is then to propose an alternative to this method in order to approach the Q.E.S. equations 
with a known number (we will concentrate essentially on three) of analytic solutions in a general way. To do so we 
will make use of the general algorithm developed in |^ . 

The paper is organized as follows: In Section wc explain the formalism of our method. It is illustrated through 



examples in Section [II. In Section [V, we compare our method with the one of Kuliy and Tkachuk. Finally, we draw 
our conclusions in Section |^. 

II. THE GENERAL APPROACH 

We concentrate on the one-dimensional and time-independent Schrodinger equation 

iJNix) = ENlpNix). (1) 

The eigenfunctions %1^n{x) as explained in are written as 

N 

i,^{x) = g{x) [f{xt^ ^™ ^ [M^)]" (2) 

where the function g{x) is a weight factor, [/(a;)]^" is introduced due to eventual singularities in V{x), and the 

quantities (real or complex) Cm ^ are expansion coefficients on the basis h{x). Actually this basis plays the prominent 
role in our approach in the sense that each of the functions entering our formalism can be expanded in this basis 
according to 

M 



g'{x) = -g{x)Y,gnKx)\\ (3) 

M 

/(^) = E/°[M^)]'' (4) 

M 

/'W = (5) 

M 

h'{x) = Y.^}[h{x)]\ (6) 
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the prime standing for the derivative with respect to x. The upper-right index for the coefRcients g] ^ ff ^ fl ,h\ 
indicates the derivative order of the corresponding expanded function. The upper summation index M is in general 
an arbitrary positive integer, big enough to be the highest value of all upper summation indexes in Eqs. the 
added expansion coefRcients being vanishing when necessary. Finally the upper summation index N is also a positive 
integer left arbitrary until we consider specific examples in the next Section. 
In 1^ , we have limited ourselves to specific potentials of the form 

V{x)= £ Vu[f{x)tY.viAh{x)]' (7) 

the coefficients Vk and vi.k being fixed according to the interaction we want to deal with while I runs in the space of 
integers. The positive integers Vmin and Vmax are also determined once the interaction is given. The expansion (|^) 
has been proposed to include all known Q.E.S. potentials. However, in the present paper, we do not have to ask that 
the potential obeys a given expansion due to the fact that the number of eigenfunctions is fixed at the start. (We will 
nevertheless see in the next Section that the determined potential will be of the type of Eq. (^).) 

More precisely, if we consider three eigenstates as stated in the Introduction, we have two possibilities: Either these 
eigenfunctions are of the form given in Eq. (||) with N = L,L + 1,L + 2, or they coincide with Eq. (||) with the same 
fixed N = L, but different values of ^ in order to have linearly independent polynomials in the basis h{x). The 
analysis of already known Q.E.S. potentials reveals that the Kuliy-Tkachuk approach is included in the first possibility 
while the Lie algebraic one can be found in the second possibility. 

In both cases, we have the first eigenstate equal to 

L 

M^) = 9i^) ifi^t' E IM^)]" (8) 
and, for simplicity, associated to the energy 

El = 0. (9) 

Such a requirement determines the potential V{x) up to the presence of parameters which will be fixed asking for the 
existence of the other solutions. More precisely, by inserting Eqs. ^ and (||) in Eq. (|I]) and using Eqs. (||)-(^), the 
potential is 

v{x) = [hixr+^^-y [fJ^gUA'^ 

+ fkf^hlh\cfh{s ~l + m) 1. (10) 



In this expression and the following ones, we have, for simplicity, omitted all the summation indexes but it is understood 
that there is a summation on each repeated index (except L). 

Now the potential of Eq. (|l^) has to be constrained so that there exist other solutions to it. For the first possibility 
mentioned above, i.e., the one including in particular the Kuliy-Tkachuk potentials, we just have to require that the 
functions (||) with = L -I- 1 and N = L + 2 are solutions of Eq. ([I]) with V{x) fixed according to Eq. (pX|). This 
gives rise to (two) systems of algebraic equations (instead of differential ones). This is the main advantage, besides its 
generality, of our approach based on the development of all functions in the /i(2;)-basis. These systems are specified 
by vanishing the different coefficients of linearly independent terms in 

-fkl^hlhlA'A^Mt -l + m) [h{x)r'+-+-+^^'-' + 2tfj^glhlci'^c^^ [/.(:.)] 

- 2X^f°flthlci^^ci''^ [Mx)]^-+'+'"+^+*-^ + (A^ - xJf.flhllci^H''^ [hixt^'-^-^^^'-' 



- 2f,glfWA''^ [h{x)r'+"'+^+' + (Xl + X^- l).flflci^^c[^^ [h{x)r'^'+' 

+ 2XUlllhlscf^cf^ [h{x)r'+^-^^+'-' - 2fjUlshlcf^c['^^ [Mx)]^-+'+™+"+^+*-^ 

+ fkf?hlhls{s - 1 + m)c(^)ci^) [Mx)]'=+'+'"+"+^+*-^ - E^f,fUf)c[''\h{x)f+'+^+' = 0, (11) 
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where N = L + 1, L + 2. It is clear that if we want to consider the case of two only eigenstates, we have to solve only 
one system, the one corresponding to iV = L + 1. We can also extend this approach to the case of more than three 
solutions. The potential of Eq. (|o|) is then more and more constrained until it is impossible to fulfill the constraints. 

In what concerns the second possibility, i.e., the one associated to the Lie algebraic approach, we have to ask for 
the function in Eq. (|^) to be a solution of Eq. with N 
The system of constraints is now: 



L, together with other coefficients than Cm\ say Cm\ 



-f.hlhlci^^dl'^^t - s)it -l + m + s) [hix]] 



/c+m+n+s+i — 2 



0. 



(12) 



The number of independent sets of coefficients c[^^ satisfying Eq. ( p^ ) gives the number of eigenstates, besides ipLix) 
given in Eq. (||). 

III. EXAMPLES 

Let us now illustrate our approach through specific examples. The point is simply to choose the function f{x) in 
an arbitrary way and then to adjust the important basis function h{x) so that Eqs. are satisfied. Our recent 

analysis of Q.E.S. equations in ||] has revealed that the choice 



f{x) = l + : 



(13) 



covers many Q.E.S. interactions. This, as well as f'{x) 
Eqs. (0)-(B|)i get the non vanishing coefficients 



2x, clearly suggests the basis function h(x) = x. From 



fS = lJ^ = lJl=2,hl = l. 



(14) 



For the first possibility, corresponding to the constraints in Eq. (|ll[), let us fix L = 1 and choose, without loss of 
generality, c^^-* = 1. The potential of Eq. ( p^ ) then reads 



V{x) = 9l,glx"^+- 



WmX ^ i 



1 + 



4Aia; 



„1 „m 



(1+^2)(,(1)+^) 



(1) 



(l + a;2)2 



2 i'm'^ 



(15) 



A rapid look at the first constraint, corresponding to = 2 in Eq. (|1 Ij) , shows that the running index m has in fact 
to be restricted to the values to = and m = 1. This first constraint is then a system of eight (nonlinear) equations 

for the eight unknowns g}),g\, c'^\ c'^K c'i \ Ai, A2, E2 if, once again, we take advantage of the fact that the eigenstates 

(2) 

are fixed up to a normalization constant, so that can be chosen equal to 1. The second constraint, corresponding 
to = 3, is a system of nine (nonlinear) equations for the nine unknowns ^q, 17]^, Cq^\ Cq'^\ c^'^\ Cj'^'', Ai, A3, i?3 with, 

(3) 

for the same reason, c\ =1. Thus there is a total of seventeen equations for thirteen unknowns. Nevertheless, this 
system is rather easy to handle and furnishes two sets of solutions. The first set is the one corresponding to the 
well-known harmonic oscillator potential while the second one is 



lo = 0, 31 = Vc, Ai = A2 



Vc- A3 = ^ - Vc, 



E2 - 2^c, E3 ^ -8c + 12y/^, c[, 



(1) 



-l,ci- 



(2) 



.(3) 



(3) 



(16) 



with c being 



■ arctan 



(17) 
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i.e., c = 1.119639. The associated potential is, from Eq. (10), 

V{x) = cx^ - 4c - Vc + — ^ — 

^ ^ 1 + a;2 (1 + 2:2)2 

and is, to our knowledge, a new Q.E.S. potential. It admits three analytic (and physical) solutions given by 



■01(3:) = xil + x^) 



exp 



= (l-.i:2)(l + a;^) 



exp 



Ei=0, 

2 - 



/CX 



E2 - 2^, 



l-l2--c\x' 



(l + x^) 



2n-v^4 



exp 



/CX 



Es = -8c+ 12Vc. 



(18) 

(19) 
(20) 
(21) 



We remark that the same analysis, but with i = 0, would have led to the potential given in @. Thus, our method 
includes the one described in that reference. Moreover, as we will show in the next Section, it is simpler (see the 
function U{x) of Eq. ( ^ which would be needed to find the potential E g. (H^ ) and the eigenfunctions Eqs. (pj|)-(^). 
For the second possibility, corresponding to the constraints in Eq. (|l2|), let us fix L = 2, which is the minimal 

value for having at most three eigenstates in this case, together with c]^^' = 1 and 



1. We have not found a 



new potential with three eigenstates. Only known ones like the harmonic oscillator or its supersymmetric partner ||] 
emerge. We have however found a new potential with two eigenstates. It is given by 



V{x)^ 



and the two eigenstates and related eigenvalues are 



1 



2ix 



36 6 3(x + i)2 



Mx) = i3 + 2ix + x'){l + x')-'exp{-—],E2 = 0, 



= i-l + 2ix + x^)cxp{-—],E2^- 



(22) 

(23) 
(24) 



Thus, we have here an example of a complex but PT-symmetric Q Q.E.S. potential, with real eigenvalues. This 
potential is new compared to the ones already studied through other developments |p^ . 

To end this Section, we mention that the potentials Eq. (^Sj) and Eq. ( p2| ) could be written as in Eq. (|^), implying 
that, given the potentials, we could have found the eigenfunctions using the algorithm we have described in |q]. 



IV. COMPARISON WITH THE KULIY-TKACHUK APPROACH 



For completeness, let us briefly recall what the method of Kuliy-Tkachuk is. It is based on a recursive construction 
of supersymmetric partners. Indeed, the initial potential V{x) is required to be written as 



V{x) = WL{xf ~-Wi{x). 
If such a Wl{x) does exist, then the ground state and its energy are known 



iPl{x) = exp 



WL{x)dx 



Er = 0. 



Provided one can find two other functions Wl+i{x) and Wl+2{x) such that the Ricatti equations 

{x) + Wl {x) +Eu^ Wl^^ {x)^W;,+,{x)+Ek+i, k^L,L+l 
are satisfied, two other solutions (respectively associated with -El+i and El+2) are available. They are given by 



1pL + l(x) = 
i^L+2{x) = 



-^ + Wl{x) 
ax 



exp 



^— + Wl+i{x) 
ax 



WL+i{x)dx 
exp 



WL+2{x)dx 



(25) 

(26) 

(27) 

(28) 
(29) 
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The main point of the Kuhy-Tkachuk approach is to reduce Eqs. <\27i} to only one equation 



{L+l), 



(x) 



Wi^\x) 



d 

dx 



{El+2 - El+i) W^f Hx) + EL+iWY'^'>ix) = 



Wl^\x)wl^+^\x) 

(L+l), 



of the two functions 



wi^\x) = Wk{x) + Wk+i{x), k = L,L + l. 



FinaUy, Eq. (M) can be considered as an algebraic one in terms of w\^\x), for instance, if the function 



U{x) = w\_^\x)W, 



(L+l) 



(x) 



(30) 



(31) 



(32) 



is fixed arbitrarily at the start. The whole approach thus relies on an adequate choice of this function U{x). Note 
that the authors of Ref. Q have restrained themselves to the context of the functions U{x) such that 



U{x) <0 if X e [min(a;o, io), max(xo, So)]) 
U (x) > otherwise 

where xq and xq are the two zeroes of U (x) such that 

?7'(min(xo, xq)) < 0, C/'(max(xo, xq)) > 0. 



(33) 
(34) 



(35) 



In 1^, these conditions were fixed according to the requirement that none of the functions Wi{x), I — L,L + l,L + 2 
possesses singularities. 

Now we concentrate on the determination of this U{x) related to the examples developed in the previous Section. 
In what concerns the potential Eq. ( p^ and its eigenstates (p^-(pTl), a comparison with Eqs. ( p6| ) and (28)-([29|) 
gives rise to the identifications 



(36) 
(37) 
(38) 



Wl{x) = 


Vex 


WL+lix) = 


^/cx 


Wl+2{x) = 


\fcx 



(2V^ 


- l)x 


1 








(1 + 


.x^) 


X 










+ l)x 


1 








(1 + 


x^) 


+ -, 
X 










+ l)x 


1 


2x [(14c- 


3V^- 15 


) + (68c+ 12Vc- 90)x2] 


(1 + 


X2) 


X 


(34c + 6Vc 


- 45)x4 + 


(14c- 3Vc- 15)x2 - 3Vc 



It is already evident that such functions will lead to a Q.E.S. potential excluded in the Kuliy-Tkachuk approach due 
to the singularity x = in Eqs. (^6|)-(^8|). Indeed the corresponding U{x) can be obtained through Eqs. (31)-(|3^) 
and is 



V{x) = 2^fcx^{x'^ - 1) 



(-192c - 39V^ + 255)x2 + 2(-17c + 3^^ + 15) 
(34c + 6^/c - 45)x4 + (14c - 3Vc - 15)x2 - i^c 



(39) 



By opposition to the U(x) permitted in |^], this one admits one second-order zero point (x = 0) and two ordinary 
zero points (x — ±1). Moreover U{x) has no pole and is such that (compare with Eqs. (^3|)-(^) 



Uix) < if X e [-1,1], 

C/(x) >0 if X e] - oo,-l[U ]l,+oo[ 

[/'(O) = 0, c/'(-i) < 0, C/'(l) > 0. 



(40) 
(41) 
(42) 



Let us also stress the heaviness of the U(x) given in Eq. (|39|): The least one can say is that it is not natural to think at 
the start of such a function while the choice of our function /(x) = 1 + x^ is on the other hand quite straightforward. 
This is the main advantage of our approach: using L = as a starting point in Eq. (|^), we could have obtained the 
Kuliy-Tkachuk solution of Ref. [^ doing the same type of easy calculations that we have performed here. We would 
have of course obtained the same results as has been used in that reference, i.e.: 



V{x) 



3 2 



l(6-7V3)- ^(f + f) + ^H + 2V3) 

2 V / l-f x2 (1-1- x2)2 



(43) 
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M^) = (1 + exp (~^xA , i?o = 0, (44) 



4 

^-2(2;) = x{l + x2)(v^-i)/2 exp (^-^^^^ , i?i - 6 - 3V3, (45) 
^-3(2;) = (l-x2)(l + a;2)(^-i)/2exp j-^x^j , £^2 -6-2%/3. (46) 



On the contrary, it is doubtful that the function U{x) of Eq. ( |39[ ) could have been witnessed from the start, while, in 
our approach, it just consists in starting with L = 1 in Eq. (^). It is also interesting to remark that the potential V{x) 
of @, and our potential of Eq. ( [l8| ) have the same type of x-dependence: V{x) — a + bx'^ + c/(l + x'^) +d/{l + x^)^, 
and admit three solutions. This shape of the potential is very rich, see e.g In fact it is even possible (apart from 
the harmonic oscillator c = d = 0) to find an exactly solvable model, using a = Q,b — 1/4, c = 4, d = —8. One may 
wonder why it is not just possible to study the Schrodinger equation with a potential of this general shape. We have 
shown in Ref. j|] that the constraints are extremely difficult to handle, especially when the number of solutions is 
growing. 

For the potential of Eq. (p2|), the function U{x) simply reduces to w['"\x) as defined in Eq. (|3l]) due to the fact 



that this potential has two eigenvalues only. We respectively have 



evidently leading to 



^ 6 X + i ^ X + 3i x + Ai ^'^ 



!/(.).<)(-) = I- (49) 



So once again singularities (in the complex plane) arise and the resulting U (x) could not be found following the 
method of Ref. §]. 

Remarks concerning the choice ofU{x) 

We have already seen that the choice of U (x) is a non trivial task. We want to show here that the choice is also 
sensitive to small changes of parameters, and that care must be done to extract meaningful solutions, i.e., even if U{x) 
allows to generate a well-behaved potential V{x), the eigenfunctions obtained using the formalism are not physically 
suitable. We show this starting with the simple choice U{x) = cx^, with c a constant to be determined. To avoid a 
0(1/2;^) singularity in the potential, one is forced to adjust the constant c to c = Ei {E2 — Ei), with Eq — 0. Choosing 
El = 1,E2 — 2, so as to have an equally spaced spectrum, two potentials emerge: either V{x) = (x^ + 10)/4, or 
V{x) = (x^ — 2)/4, which are both harmonic oscillator potentials. Because there is only a translation between the 
two, one may believe that both are admissible. This is however not the case using the general method of Ref. 
Using the different functions Wl,L = 0, 1,2 which enter Eqs. (Pq)-(p2), the case which would lead to V{x) implies 



unbounded solutions Eq. ( |26[) , Eq. ( pSD and Eq. ( [291) at infinity. Only the case which lead to the potential V{x) has 
WL-functions which ensure that the eigenfunctions are physical. Moreover, changing minimally the problem leads 
often to the impossibility to solve it. For example, taking U{x) — 3x^, corresponding with non equally spaced energy 
levels Eq = 0, £^1 = 3, i?2 = 4, leads to well-defined potentials, because Wq{x) behaves adequately. It is however 
impossible to find physically acceptable wave- functions due to singularities and bad behavior at infinity of Wi {x) and 

W2{X). 

Another remark concerns the class of functions given in Ref. which is too restrictive. For example, the well-known 
Q sextic oscillator 

Vb(x) =x^-llx2 + 8 (50) 

admitting the three solutions 

iPoix) = (l + 4x2 + 2x'')exp(-ixM , £;o = (51) 
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Mx) = (2x4 - 3) exp ( -^a;4 j , Ei = 8 (52) 

„2 , o„4x „„„ I 1 4 



corresponds to 



i/,2(x) = (l-4x^ + 2x^)exp ( --x"* ) , £2 = 16. (53) 



This U{x) does not enter the Kuhy and Tkachuk formahsm. 

We have also tried different other choices for U{x), which do not belong to the class of Kuliy and Tkachuk. We 
have investigated the two simplest deviations of the quadratic choice: U{x) — c{x^ — a^) and U{x) — cx^ix^ — a^), 
with both choices > or < 0. The result is that it is possible to find solutions for the first case. These solutions 
are however not physical (they diverge at infinity). In the second case, it is possible to find physically acceptable 
solutions. However, it is not possible to express them analytically because the functions Wl,Wl+i and Wl+2 which 



enter Eq. (|26|) and Eqs. (E8h-(29) are not analytically integrable. 



V. CONCLUSIONS 

In this paper, we have presented a powerful alternative to the method of Kuliy and Tkachuk of Ref. We 
have shown that the function U{x) is, even in a case which differs minimally from the one treated in that reference, 
complicated and heavy to handle: simple U{x) do not lead to Q.E.S. potential; other U{x) may give a good potential, 
but not physically relevant solutions. In contradistinction, our method seems straightforward to apply and, by 
construction, leads to physical solutions. 
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